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1. Introduction 

Chiral perturbation theory (chPT) ^ ||] is a widely used tool in many phenomenological 
applications and also helpful to guide an extrapolation to lighter quark masses in lattice-QCD 
simulations. Here we will report on a determination of the NLO low-energy constants (LECs) Ij, 
and /~4 which appear in the light quark mass dependence of the pseudo-scalar meson masses and 
decay constants in SU(2) chPT. 

We analyze configurations generated by the Wuppertal-Budapest Collaboration [Q, ||, |6[ [7], ^ 
using the Symanzik glue and 2-fold stout-smeared staggered fermion action for a 2+ 1 flavor QCD- 
simulation. The mass of the single flavor has been kept at the value of the physical strange quark 
mass, whereas the two degenerate lighter quark masses have been varied such that light meson 
masses in the range of 135 to 440 MeV were simulated. The simulations were performed at five 
different gauge couplings /?, resulting in lattice scales between 0.7 and 2.0 GeV (see next section 
for details on how the scale has been determined). Table [j] summarizes some of the parameters of 
the simulations. 

The 2-fold stout-smeared version of the staggered quark action has been proven to be advan- 
tageous [^] in reducing the inevitable taste-breaking of staggered fermion formulations. Therefore, 
in this work we only consider the pseudo-scalar mesons with taste matrix 75 when measuring me- 
son masses or decay constants. Details of the computation of these quantities will be reported in a 
forthcoming publication. 

2. Scale setting and physical quark masses 

To set the scale at each simulated gauge coupling j8 and identify the physical point, i.e. the 
average up/down quark mass m^ hys = (m u + md)/2 corresponding to a pion in the isospin limit 
with an estimated mass of M n = 134.8 MeV we use a two-step procedure. First, we ex- 
trapolate the ratio (aMu) 2 / (afu) 2 of the squared meson masses and decay constants to M 2 /f 2 = 
(134.8MeV/130.41MeV) 2 = 1.06846, where we also used the PDG- value = 130.41MeV [jfij|. 
In that way arrtf hys is obtained. In the second step, we extrapolate afu to this quark mass value and 
obtain the lattice scale with the help of the PDG-value for f n . For the extrapolation we used two 
different ansatze: a quadratic and a rational (linear in numerator and denominator) fit form. An 
example of these extrapolations is shown for the ensembles at j3 = 3.85 in Fig. [l[ There, like for 
all other j3-values as well, the heaviest quark mass point has been excluded, resulting in a fit range 
of approx. ami/am^ hys < 8.0 (corresponding to Mu < 390 MeV). We stress that here, like in the 
chiral fits to be discussed below, the data has been corrected for finite volume effects beforehand, 
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1.0, 4.0, 6.0, 8.0, 10.0 


40 3 x 64 


- 16 3 x 32 


3.85 


2.04 


1.0, 1.4, 2.0, 4.0, 6.0,8.0, 10.0 


48 3 x 64 


- 24 3 x 48 



Table 1: Simulated lattice ensembles: gauge coupling j3, lattice spacing 1/a, simulated quark masses m\, 
and range of lattice sizes. 
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by means of using the two- and three-loop resummed formulae of [11] for the pion decay constants 
and masses, respectively. Our spatial lattice volumes L 3 are in the range (4.7 fm) 3 - (6.8 fm) 3 with 
a minimal MuL f» 3.2, ensuring that the finite volume collections within our fit ranges are at most 
at the order of 1 per cent for the decay constants and even less for the meson masses. 

By fixing I /a and am^ hys in the way described above, the meson masses and decay constants 
show no discretization effects at all directly at the physical point and we can assume those effects 
to be small (since of higher order in the quark masses and/or lattice spacing) in the vicinity of the 
physical point, i.e. in the mass range covered by our fits. Such discretization effects, of course, are 
present in other observables, which are not considered in this work. 



3. Fits to NLO SU(2) chPT 



The quark mass dependence of the finite-volume corrected data for the meson masses and 
decay constants is fitted simultaneously at different j3 -values using the NLO-SU(2) chPT formulae 
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where we made use of the already determined 1 /a and am^ hys to scale the quark masses and the 
meson masses and decay constants measured in lattice units. This fit has four free parameters: two 
NLO low-energy scales A3, A4, the decay constant in the SU(2) chiral limit / and the renormaliza- 
tion scheme-independent combination (2Bm^ hys ) of the LO low-energy constant B and the physical 
quark mass m^ hys . 

We would like to point out that the chiral fit formulae do not include any taste breaking effects, 
i.e., we did not use staggered chPT. This seems justified to us, since we are only considering 75- 
taste mesons as mentioned above and use these to define our scaling trajectory at the physical 
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point. In other words, since the meson mass and decay constant at the physical point were used to 
set the quark masses and lattice scales, no discretization or taste breaking effects are present in the 
chPT formulae for M 2 , and fu as discussed above. Furthermore, taste breaking effects are reduced 
anyway by the choice of the fermion action as mentioned above. 

The top panels of Fig. || show the combined fits including the data at all lattice spacings and 
for meson masses in the range of 135 to 390 MeV. (Here and for all following plots we mark data 
points included in the fit by circles, while those not included in the fit are marked by diamonds.) 
As one can already see by eye, the description of the data by the fit is not satisfactory, resulting in 
a # 2 /d.o.f. « 4.3. As expected, the fit quality measured, e.g., by # 2 /d.o.f. improves continuously 
when reducing the upper bound of the meson mass range. The middle panel of Fig. |2| shows the 
fit to all meson masses in the range 135 MeV < Mu < 275MeV giving an acceptable # 2 /d.o.f. py 
1 .0. A similar improvement can be achieved by excluding the two coarsest lattice ensembles from 
the fit, i.e., limiting 1/a > 1.3GeV. The bottom panel of Fig. |2| shows an example of such a fit 
with 135 MeV < Mu < 340 MeV resulting in # 2 /d.o.f. «1.7 (this number has to be compared 
to # 2 /d.o.f. « 2.6 for the same mass range and using all j8). Applying both kinds of cuts, i.e. 
135MeV < Mu < 275MeV and 1/a > 1.3GeV, eventually gives a £ 2 /d.o.f. ps 0.8. 

Here we are mainly interested in the SU(2) low-energy constants Z3 and I4 which are related 
to the low-energy scales A3, A4, respectively. Therefore, in Fig. || we show the fitted values for 
these parameters obtained with different fit ranges. We also display the ratio f K /f as obtained 
from the various fits. Whereas for % (left panel), if at all, one could identify a shift in the result 
depending on whether or not the two coarsest lattices ensemble are excluded, for 1 4 and f n j f one 
observes a clear dependency on the fitted mass range, while the influence of excluding coarser 
lattice ensembles seems to have only a marginal effect. Eventually, we quote as our result for the 
low-energy constants the central value and statistical error obtained from the fit range 135 MeV < 
Mu < 275 MeV, 1/a > 1.3GeV and take the variation with respect to that value from other fits 
including the nearly physical points (data marked by asterisks in Fig. ||) as our estimate for the 
systematic error, so that we obtain: 

h = 2.90(1 l)stat(17) S yst, h = 4.04(04) stat (13) syst , f n /f = 1.0627(07) stat (24) syst . (3.4) 

Since often lattice data from meson masses larger than the physical M n are extrapolated to the 
physical point using SU(2) chPT, we also investigated fit ranges excluding the physical point. In 
Fig. [I] the fits for the meson decay constant are shown for 230 MeV < Mu < 340 MeV (left panel) 
and 230MeV < Mu < 390 MeV (right panel). As one can see in the close-up view of the region 
near the physical point, the value for f„ extrapolated from such a fit is below (/£ xtr ' ~ 128(1) MeV) 
the values simulated near the physical point. As one can see from Fig. [3J, also I4 and f K j ' f are 
significantly changing, once the nearly physical points are excluded from the fit range. 

4. Fits to NNLO SU(2) chPT 

Extending the SU(2) chPT fit formulae for the meson masses and decay constants to NNLO 
(e.g. cf. |]|]), in our set-up three new fit parameters have to be added: a combination of the NLO 
low-energy constants l\, l?. In = (Jh +8/2)/15 and two parameters for NNLO-LECs k m , kf. 
Again fitting our data for the meson masses and decay constants at various j8 simultaneously now 
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Figure 2: Combined NLO SU(2) chPT fits with various fit ranges. The left panels show the decay constants 
///, the right panels the squared meson masses Mf, divided by the quark mass ratio ami / a>rf hy * . The fit 
ranges are: top: all j8 , 1 35 MeV < M„ < 390Me V, middle: all j3 , 1 35 MeV < M u < 275 MeV, bottom: only 
1/a > 1.35GeV, 135MeV <M„ < 340MeV. 



using the NNLO fit formulae without any constraints on the fit parameters (7 in total) leads to an 
unnatural order of the NLO- compared to the NNLO-contribution as can be seen from the left panel 
of Fig. There the black line denotes the full fit up to NNLO and the red line only the contribution 
up to NLO, the large difference between the two being the NNLO-contribution. The situation can 
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Figure 3: LECs obtained from NLO SU(2) chPT fits with different fit ranges: left panel: 1$, middle panel: 
U, right panel: f%j f. Blue points denote fits where l/a> 1.35 GeV, red points fits where all j3 are included. 
Fits including the nearly physical points are marked by an asterisk. The solid, dashed and dashed-dotted 
lines display the central value, statistical and combined (stat. and syst.) error, resp., of our quoted results. 





Figure 4: Two examples for NLO SU(2) chPT fits excluding the nearly physical points (only /// is 
shown here). Left panel: 230MeV < M u < 340MeV, right panel: 230MeV < M u < 390MeV; both 
l/a> 1.35 GeV. 

be improved by using priors for some of the fit parameters, e.g., using a phenomenological estimate 
for l\2 = 2.1 ± 0.3 as can be obtained from values quoted for fi, I2 in [ 12]. A fit using such a prior 



is shown in the right panel of Fig. BL which describes the data well and has a reasonable ordering 
of the NLO- compared to NNLO-contribution. Still we refrain from using NNLO-chPT as long as 
we do not have enough data in the light quark mass region to constrain such fits without having to 
rely on additional input used for priors on the fit parameters. But it is reassuring to us, that by using 
such priors, a NNLO-fit results in NLO-LECs comparable to those found in our NLO-fits. 



5. Conclusions 

From our NLO SU(2) chPT fits to meson masses and decay constants measured on staggered 



2+1 flavor lattice simulations of QCD, we quote the following set of LECs (see Eq. (p.4|)) as our 
preliminary result: 



h = 2.90 ±0.20, U = 4.04 ±0.14, f % jf = 1.0627 ±0.0025 . 
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Figure 5: Two examples for NNLO SU(2) chPT fits (only /// is shown here). Left panel: without priors, 
right panel: using prior l l2 = 2.1 ±0.3. Both use all j5 and 135MeV < M u < 340 MeV. 

These values are in good agreement with other recent lattice determinations of LECs, for 
example the FLAG-report [[J quotes h = 3.2 ±0.8 and f K /f = 1.073(15) as lattice averages, 
while due to some tension in the results no value for I4 is quoted at the moment. Our findings 



also agree well with the phenomenological estimates Z3 = 2.9 ±2.4 and I4 = 4.4 ±0.2 [|12j and 
Ulf =1.0719 ±0.0052 flT|g]. 

For a forthcoming publication we hope to have additional data points available at light quark 
masses corresponding to meson masses between 135 and 275 MeV. More details about our chiral 
fits will be reported there as well. 

The speaker acknowledges support from the DFG SFB/TR 55 and the EU grant PITN-GA- 
2009-238353 (ITN STRONGnet). 
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